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Abstract 

We explore the consequences of Bose-Einstein condensation on two-scalar-exchange 
mediated forces among bodies that sit in a boson gas. We find that below the 
condensation temperature the range of the forces becomes infinite while it is finite 
at temperatures above condensation. 



Van der Waals type forces where two photons are being exchanged |1[ or the extremely 
feeble forces generated by 2-neutrino exchange |3| provide examples of forces among two 



static bodies in a vacuum produced by the exchange of two quanta in the t-channel. Spin 
independent interactions arising from double (pseudo) scalar exchange 0, [| || (such as axions 
and/or more bizarre specimens of modern completions of the Standard Model(SM)) provide 
further examples for these so-called dispersion forces 0. When the objects that feel such 
forces are placed in a heat bath at a temperature T, the forces get modified. Indeed, in the 
case of molecules in the relic photon background, the long-range Casimir- Polder forces among 
them are strongly affected for distances much larger than T _1 || and, for the 2-neutrino forces, 
the cosmic neutrino background completely screens off the interaction at large distances 
(again, large meaning much larger than T" 1 ). 

In the present paper we shall deal with a gas of scalar bosons carrying an abelian charge and 
non-zero chemical potential. As mentioned before, their double exchange between fermions 
has been studied in vacuum. The case of a non-charged scalar bath in a classical Boltzmann 
distribution was briefly discussed in ||. However we are not aware of discussions on the effects 
resulting from placing the interacting system in a charged scalar heat bath displaying genuine 
quantum statistical effects such as Bose-Einstein condensation. Because in the previously 
reported instances, interesting effects did result, we think it is worthwhile to raise this issue 
here. Admittedly, light scalar bosons have a much different status than photons and neutrinos 
and their nature is entirely speculative. Nonetheless, in almost any extension of the Standard 
Model, scalars are present and some have been furthermore suggested as candidates for dark 
matter so that might be part of the cosmic relic background. Not to speak about axions 
needed to solve the serious CP problem of QCD and sometimes advertised also as galactic 
halo dark matter or even as dark matter on a cosmic scale. 

Since we do not pretend to resemble any particular model of matter scalar interactions we 
will construct a very simple toy model that avoids unessential complications and that might 
mimic more "realistic" interactions of matter and light scalar fields. Consider the interaction 
lagrangean, 



$ is a heavy scalar field of mass M and (p is a light scalar field of mass m. Let us now put 
two such heavy particles $ in vacuum at a distance r. Their lowest order interaction is given 
by the Feynman amplitude in figure 1. The potential is obtained from the NR limit of this 
amplitude via Fourier transformation. That is, 
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Figure 1: Feynman diagram corresponding to the interaction which gives rise to the long-range 
force. 

where ^ 

T ^ = 9 I (2vr) 4 k 2 - m 2 + ie (q - k) 2 - m 2 + ie ^ 
Taken at face value the previous integral diverges. It has to be regulated and the piece 
which leads to a long-range interaction extracted. We are left with a function of logQ 2 , with 
Q = |Q|, after dropping any polynomial in Q that would only yield contact interactions |Tt] 



Introducing a convergence factor e v ® to evaluate the Fourier transform and taking the limit 
f|^0we obtain: 

for the potential. The same result can be derived using the dispersion theoretic techniques of 
Feinberg and Sucher 0. 

For small <p mass, eq. (^) becomes: 

g 2 

Kac(r) = ~128vr3r3M 2 (5) 

valid for r 1/m. Beyond this range the Bessel function gives rise to the characteristic 
Yukawa factor e~ 2mr . 

Incidentally, identical r behaviour obtains for the spin independent potential arising in 
double exchange of pseudoscalars coupled to matter fermions via Yukawa couplings || |J . 

Next we introduce our system in a heat reservoir made of an ideal relativistic (p gas at 
temperature T (T > m). We further assume that the particles in the gas carry a conserved 
quantum number (which we will refer to as " charge" ) corresponding to a quantum mechanical 
operator Q. We may use real time finite temperature field theory |TT|] to calculate the effect 



of the heat bath on the potential between the two massive particles. So, we take for the 
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T-dependent (^-propagator, 

1 

k 2 — m 2 + ie 



D F (k,T) = — 9 \ - 2m5{k 2 - m 2 ) \9(k°)n(\k°\,T) + 0(-k o )n(\k o \,T)] (6) 



where 

n(u, T) = (e {w -^ /T - l) ~* and n(uj, T) = ( e ( w +">/ T - l) ~* (7) 

are the B-E distribution functions for particles (charge +1) and antiparticles (charge -1), 
respectively. \i is the chemical potential associated to the conserved charge Q. 

The amplitude of figure 1 now generalizes to 

T(q) =9 2 J ^0^^ (k + k'- q)D F (k, T)D F (k\ T). (8) 

This amplitude generates two distinct contributions to the potential. The first one arises from 
the first piece in D F (k,T) and is the vacuum potential just derived. The other corresponds 
to the situation where one of the scalars in the double exchange process is supplied by the 
thermal bath. This effect is described by the crossed terms in the amplitude involving the 
thermal piece of one (^-propagator along with the vacuum piece of the other propagator. This 
thermal component of the Feynman amplitude can be written as 



T T {q) = -ig 2 J 



d 3 k 1 



(27r) 3 2 v / k2T 



1 1 



(n + n) (9) 



(k — q) 2 — m 2 (k + q) 2 — m 2 

where use has been made of energy-momentum conservation and of the S(k 2 — m 2 ) in D F (k, T). 
In the static limit, i.e. momentum transfer q ~ (0, Q), where matter is supposed to be at rest 
in the frame of the heat reservoir, the piece in square brackets in eq. (^|) reads: 

-2(Q 2 -4k 2 Q-k)~\ (10) 

So, finally the T-dependent amplitude to be Fourier transformed is 

r T fa g (o, Q))-*/^ Vk2 + m2 Q2 _ 4k2Q . k (» + n) (ID 

i.e. it has been reduced to an integral over the phase space of the real particles (and antipar- 
ticles) in the heat bath. 

The reservoir is thermodynamically characterised by a temperature T, a volume V, and a 
fixed charge QQ. Then, the chemical potential fi(T) is determined from the relation 

Q = Y,{n-n). (12) 

k 



1 For definiteness we take Q > 0, i.e. particles outnumber antiparticles. 
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For a Bose-Einstein gas, the sum over states in the previous formula can be converted to an 
integral like the one in eq. ( |TT| ) as long as its temperature is above a critical temperature T c . 
Below that temperature, if one makes the replacement 

d 3 k 



in eq. fll2|) , the result is less than Q fL2| . This is because below T c a large macroscopic 



fraction of the charge resides in the lowest energy state and the density of states Vk 2 /2tt 2 in 
the continuous representation of the sum over states gives a zero weight to the zero mode. 
On the contrary if the gas is above T c then the charge is thinly distributed over the states 
and no individual state is populated by a macroscopic fraction of the total charge so that 
by passing to the continuum essentially only an infinitesimally small error is done. Let us 
discuss both cases in turn. First start with the non degenerate case, i.e. when T is above 
the condensation temperature T c . In this instance, the phase space integrals in the formulas 



above correctly describe the physics of the problem. Therefore we can take eq. (|TTD and plug 
it in the expression for the potential (0) 

V T (r) = 7-7^/ dQQ sin Qr / — — -—.( n + n \ (14) 

TK ' 27r 2 r4M 2 io VV V J (2tt) 3 Vk 2 + m 2 Q 2 - 4k 2 Q • k 1 ' 1 ' 

After trivial integration over Q and over the polar angle in k-space we get 

1 



Vr(r) = <F 



g 4M 2 7r 3 r 2 



k 2 dlk 


sin 


2|k 


r 
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Vk 2 + m 2 


k 







+ 



3 (Vk 2 +m 2 +M)/T 



(15) 

This equation by itself is not sufficient to determine Vr(r) for the functions n and n contain 
the chemical potential /i(T) which has to be obtained through 



P 



f = ^/k 2 rf|k|(n-n) (16) 



which, by the way, also determines the critical temperature T c via the implicit equation [12 



p = p(T = T c ,fx = m) . (17) 

So eqs. (p~5|) and (|16D give the solution to our problem. We can use the high-temperature 
expansion of eq. ( |i~6"D derived in |T^] to obtain the chemical potential as a function of T. To 
leading order, fi(T) = m(T c /T) 2 , and we introduce it in eq. (|l"5"D to get our potential above 
the condensation temperature. The integral in eq. (|15j) can be easily split in three pieces upon 
using the following identity: 



The resulting integrals can be found in [13]. The — | piece exactly cancels the vacuum 
potential (01) so that the total contribution to the potential is finally: 



n 2 1 

irtotal y —T 1 

T ^ Tc " 64M 2 7r 2 r 2 



(19) 



fe=i 

„2,^2 i ™2/i f4\ /T^2 



with a = y (4/c7rmT c 2 /T 3 ) 2 + [m 2 (l — £ 4 ) /T 2 + 4A; 2 7r 2 ] 2 , = 4£r% 2 + m 2 (l - C)/T 2 , and 

£ = r c /r. 



Notice in the first term of eq. ( |T9| ) the typical Yukawa damping factor that cuts off the 
interaction at distances long compared to the Compton wavelength of tp. Furthermore, since 
for any k 



T v /2(a + p)> 2my/l - ^ (20) 

all modes in the second term of eq. (|T9|) are even more suppressed at large distances. This 
makes sense physically and is very convenient when trying to do a numerical evaluation. 

Below T c , a macroscopic fraction of the charge carried by particles in the reservoir piles 
up in the zero mode state (the condensate) and the integrals in eqs. ( JT5| ) and (|TB|) no longer 
correctly describe the physical situation. Indeed, eq. (|i~6D gives the density of charge in excited 



states p*, i.e. the thermal modes [12 . For a relativistic boson gas, 

p* = l -mT 2 (21) 

where we used p = m in this temperature regime since p has to be always less or equal than 
m and it monotonically increases as the temperature decreases until it reaches m at T c (and 
stays fixed in the macroscopical sense thereafter). Because the definition of the condensation 
temperature eq. (|i~7D implies in this case 

P = > (22) 



the charge density in the ground state is 

* = ' t 1 - [wj ) ■ (23) 

The Feynman amplitude eq. (ITT]) , which involves also a sum over states, should be split 
accordingly in two parts. The zero mode term on the one hand and on the other hand the 
integral over thermal modes. The zero momentum mode contributes to the amplitude 

^■(0)U = ^i(n + n)U (24) 
where the distribution function factor can be rewritten as 



since fi = m — 0(T '/ Q). 

As long as the net charge Q is a macroscopically large number many orders of magnitude 
larger than T/m, this factor coincides essentially with the condensate contribution to the 
density po- One may gain intuition on how charge is distributed among states by making a 
few numerical exercises with our formulae. By way of example, we take a fiducial volume of 
10 m 3 filled with 400 units of charge per cubic cm (i.e. numerically equal to the photons in the 
MWBR). Then, for m = 10~ 6 eV, 4 x 10 11 particles and 15 antiparticles populate the ground 
state, while 4.6 x 10 8 particles and 4.2 x 10 8 fill the excited states, at T = 0.01 T c = 3 x l0~ 4 eV. 
Clearly, the statement following eq. (23) is correct and it allows us to calculate T{Q)\k=o in 
terms of the fixed quantity p (eq. (|22"D ): 

f 
3Q 2 

The Fourier transform of this equation gives the contribution of the condensate to the poten- 
tial. It is: 

^o(r) = --^-(T c 2 -T 2 ). (27) 



^)| fc =o = ^(^-T 2 ). (26) 



48M 2 vr r 

The thermal contribution (i.e. from the excited states) is just eq. (|i~5|) with the chemical 
potential held fixed at the constant value /i = m. The same expansion, eq. (|i~8D , and the same 
integration techniques can be used now to generate the thermal potential below the condensa- 
tion temperature. Again, we find a piece that exactly cancels the vacuum contribution to the 
interaction, and what finally turns out to be the total potential, i.e. the thermal contribution, 
plus the condensate contribution (|2"T|), plus the vacuum contribution (||) is: 

MT {~ T2) + V f 1 + 2 £ e—v^) cos (2rT^ 



n 2 

-\rtotal y 

V T<T C - ^M^2 



(28) 

where 7 = kn^Jk 2 n 2 + m 2 /T 2 and r\ = k 2 ii 2 . 

It is easy to check that the potential ( ^8l) for T < T c and the potential (]i~9|) for T > T c 
coincide for T — T c (£ = 1). 

Inspection of these results leads us immediately to realize an important consequence of 
Bose-Einstein condensation. All terms in the infinite sum in eq. ( 25[) decay faster than e~ 2mr 
for T > m. Therefore, for distances much larger than the Compton wavelength of tp, i.e. 
r ^> m" 1 and hence r ^> T" 1 , V^<y is given by: 

9 2 T 2 



Vn 



total ^ y c 



rp2 rp 



(29) 



-7^2 ' r rr2' 

C c 

Namely, at low temperature (i.e. below T c ) the force, that was finite-ranged at high tem- 
perature (i.e. above T c ), becomes infinite-ranged. This comes about because the medium 
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absorbs and restores 3-momentum in the scattering process so that the 4-momentum squared 
of the other (p quantum exchanged in the t-channel can reach the mass-shell in the physical 
region of the scattering process. For the k = mode in the bath, in particular, the prop- 
agator of the second particle (see eq. (|TUD) takes the form of the Coulomb propagator, and 
becomes singular at the edge of the Q-integration region, exactly as in the Coulomb case. 
For T > T c this infinite wavelength mode has zero measure, and it does not contribute to the 
potential. However, in the condensed phase, the infinite-range potential arises as a collective 
phenomenon essentially because all charge piles up in the ground state. 

What we would like to do now is to show graphically the transition of the potential as 
we vary the temperature from T > T c to T < T c . In the first temperature regime, i.e. above 
condensation, we compute eq. fll9|) numerically and below T c we evaluate eq. fl28|) where again 
numerical methods are used. Figure 2 displays our results. 
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Figure 2: Total potential, divided by —g 2 /128M 2 ir 3 r 3 , for a relativistic hose gas made of 
particles of m ~ 10~ 6 eV and having a density p ~ 40 cm -3 . Behaviour above and below T c is 
shown. Note the Yukawa exponential damping for T > T c starting atr~l/m^rT c r^ 10 3 . 

Let us briefly summarize our findings. Light scalars are basic ingredients of many comple- 
tions of the Standard model. They may carry a new conserved quantum number. If ordinary 
matter is neutral with respect to this new charge, then these scalars should couple to ordinary 
matter pairwise. But double (pseudo) scalar exchange generates long range spin independent 
forces among bulk matter exactly as 2-neutrino exchange and 2-photon exchange (van der 
Waals forces) do. All these dispersion forces are modified when matter is introduced in a 
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heat bath. The MWBR effect on van der Waals type potentials and the relic neutrino cos- 
mic background effect on neutrino forces have been recently surveyed. The present paper 
presents an investigation of the effects of a relativistic ideal Bose gas on potentials generated 
by 2-scalar exchange. An example for such a heat bath could be provided by hot dark matter 
(i.e relativistic at decoupling) made of hypothetic relic light scalars. For this purpose we use 
a very simple model for matter-scalar interactions for we do not want to commit ourselves 
to any specific extension of the SM. We do not expect that the particular r-behaviour is in 
any sense realistic but the phenomena produced by Bose-Einstein condensation of the heat 
reservoir are totally independent of the form of the interaction chosen. What we find is a 
very dramatic effect: below the critical temperature, the finite-range force that we had above 
this temperature becomes an infinite-range force. The phenomenon arises as a combination of 
kinematics (3-momentum exchange of the matter system with the medium) and the collective 
effect of condensation of charge. In the particular model studied in this paper, a potential 
of the form ~ exp(—2mr)/r 2 at T > T c converts to a ~ 1/r potential at T < T c . Should 
hot relic scalars populate our Universe with a present density such that their temperature is 
below the threshold for Bose-Einstein condensation, then the effect described above would 



provide an excellent opportunity for experiments searching for forces weaker than gravity |L4 



for in this case no exponential decay with distance occurs and, furthermore, a milder power 
law fall-off with distance ensues. 

Work partially supported by the CICYT Research Project AEN98-1093. F.F. acknowl- 
edges the CIRIT for financial support. 
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